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The thermal conductivity of YBa2Cu30e.9 was measured at low temperatures in untwinned sin- 
gle crystals with concentrations of Zn impurities from to 3% of Cu. A linear term Kq/T = 
0.19 mW K -2 cm -1 is clearly resolved as T — *■ 0, and found to be virtually independent of Zn 
concentration. The existence of this residual normal fluid strongly validates the basic theory of 
transport in unconventional superconductors. Moreover, the observed universal behavior is in quan- 
titative agreement with calculations for a gap function of d-wave symmetry. 



PACS numbers: 74.25.Fy, 74.72.Bk 

The theory of quasiparticle transport in unconven- 
tional superconductors, developed over the last decade, 
has remained largely untested. A novel feature that arises 
when the superconducting gap function has nodes for cer- 
tain crystal directions is the existence of quasiparticles 
at T = 0. This residual normal fluid is a consequence of 
impurity scattering, even for low concentrations of non- 
magnetic impurities (see jflJ3], and references therein). 
Its presence, which should dominate the conduction of 
heat and charge at T << T c , has yet to be firmly estab- 
lished, and its properties have never been investigated. 
For certain pairing states, with appropriate gap topology 
and symmetry, an appealing phenomenon is predicted 
to occur: quasiparticle transport should be independent 
of scattering rate as T —* 0. This universal limit, first 
pointed out by Lee || for the case of a <i-wave gap in 
two dimensions, is the result of a compensation between 
the growth in normal fluid density with increasing im- 
purity concentration and the concomitant reduction in 
mean free path. 

In this Letter, we report the first observation of uni- 
versal transport in a superconductor. Our study of heat 
conduction in the high-T c cuprate YBa2Cu3 0g.9 provides 
a solid validation of the basic theory of transport in un- 
conventional superconductors and insight into the nature 
of impurity scattering in the cuprates. It also supports 
strongly an identification of the gap function as having 
e?-wave symmetry. 

The thermal conductivity k(T) was measured be- 
tween 0.05 and 1 K, for a current along the a axis 
of five single crystals: four untwinned crystals of 
YBa2(Cui_ a; Zn :E )306.9 and one crystal of YBa2Cu306.o- 
The latter was obtained by full deoxygenation via anneal- 
ing at 800 °C in helium gas for 64 h; it is insulating, with 
p a (100 K) = 42.7 ilia, x is the nominal concentration of 
Zn, achieved by mixing in ZnO powder at the start of the 
growth process in the atomic ratio Zn:Cu::1.5x : 1 — x, for 
x = 0, 0.006, 0.02 and 0.03. The experimental technique 



and the sample preparation are described elsewhere pL[| . 
The resistive T c is given in Table I. The uncertainty on 
the geometric factor is at most ± 10%, 10%, 20%, 5% 
and 10% for the x = ("pure"), 0.6%, 2%, 3% and de- 
oxygenated ("deox") samples, respectively. 

The a-axis resistivity is linear in temperature above 
130 K J§, and a fit to A + BT yields the values in Ta- 
ble I. Zn substitution has two effects: it reduces T c and 
it increases A. At low concentration, both effects are lin- 
ear, and dT c /dA = — 0.5 K / /xQ cm, in agreement with 
data on twinned crystals (e.g. 0). Concentrations of Zn 
from to 3% correspond to a large range of scattering 
rates, but to a modest level of pair-breaking: adding 3% 
Zn suppresses T c by only 20%. Given that the inelastic 
scattering term B is independent of x, the impurity scat- 
tering rate T — 1/(2t ) may be estimated via the residual 
resistivity po = m* /ne 2 To'. 



,(x) = (wjj/87r)[ P o(x = 0) + A{x) - A(0) 



(1) 



TABLE I. Sample characteristics for the four untwinned 
a-axis crystals of YBa2(Cui_ I Zn a: )306.9. x is the nominal 
zinc concentration and T c is the superconducting transition 
temperature. A and B are extracted from a linear fit of the 
resistivity (130 < T < 200 K). The scattering rate T p is 
estimated via the resistivity using Eq. (1) with u p = 1.3 eV 
and po(x — 0) = 1 pi} cm. 
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FIG. 1. a-axis thermal conductivity of the two YBa2Cu30 a 
crystals, one superconducting (y — 6.9; circles) and one insu- 
lating (y — 6.0; triangles). Main panel: k/T vs T ; lines are 
fits to a + bT 2 for T < 0.15 K. Inset: k/T vs T. 



where uo p — y/ 'Anne 2 / 'm* is the Drude plasma fre- 
quency and po(x = 0) is the resistivity of the pure 
crystal at T = 0. The latter is estimated via the mi- 
crowave conductivity, from which the mean free path 
is known to increase by ~ 100 in going from 100 K to 
~ 10 K in high-quality untwinned crystals 0. Since 
p a (100 K) ~ 75 /xfi cm, then p (x = 0) < 1 /LtO cm. 
With uj p = 1.3 eV (a-axis) one gets the scattering 
rates listed in Table I, in units of T c q—T c (x = 0). Note 
that 3% of Zn causes a 40-fold increase in T. 

In order to use k(T) as a probe of quasiparticle behav- 
ior, the phonon contribution must be extracted reliably. 
This can only be done by going to temperatures suffi- 
ciently low that the phonon conductivity n p h has reached 
its well-defined asymptotic T 3 dependence, given by 



1 



K ph = ( V ph) ^oT 6 



(2) 



where (3 is the phonon specific heat coefficient, (v p h) is 
a suitable average of the acoustic sound velocities, and Ao 
is the temperature-independent maximum phonon mean 
free path. In non-magnetic insulators, acoustic phonons 
are the only carriers of heat at low temperature and 
Eq. (2) is well verified, with (v ph ) = w L (2s 2 + l)/(2s 3 + l) 
in single crystals, where s = vl/vt is the ratio of longitu- 
dinal to transverse velocities [^). In high-quality crystals, 
A = 2w/y / 7t, where w is the (geometric) mean width of 
a rectangular sample M. 

The simplest way of investigating the phonon contribu- 
tion in YBa2Cu30j, is to remove all electronic carriers by 
setting y ~ 6.0. (Note that antiferromagnetic magnons 
are not expected to contribute at T < 1 K, since the 



acoustic spin-wave gap in YBa 2 Cu 3 06.i5 is ~ 100 K Jj-Oj.) 
The thermal conductivity of such an insulating crystal 
is shown in Fig. 1 (triangles). As seen from the linear 
fit, k/T = a + bT 2 below about 0.15 K, with a = 
and b = 14 mW K~ 4 cm -1 . The first question of inter- 
est is: what happens when electronic carriers are intro- 
duced? The answer is provided by the thermal conduc- 
tivity of a well-oxygenated crystal, also shown in Fig. 1 
(circles): electronic carriers contribute a definite linear 
term to k(T). Applying the same fit as before yields 
a = 0.19 mW K~ 2 cm^ 1 and b = 17 mW K~ 4 cm" 1 . 

It must be emphasized that such an analysis is sound 
only when applied to the asymptotic regime for K p h- To 
confirm that this is indeed the case for T < 150 mK in the 
insulating crystal, note that a = and the magnitude of 
the cubic term is right, i.e. it corresponds to a maximum 
mean free path Ao dictated by the mean crystal width 
w. Indeed, from Eq. (2) using f3 = 0.3 — 0.4 mJ/K 4 mole 
0,|| and (v ph ) = 4000 m/s (v L ~ 6000 m/s, v T ^ 
3700 m/s ||), A = 270 - 360 = 2w/y^ (see Ta- 
ble II). (Note that (3 and (v p h), given here for y ~ 6.9, 
could be slightly different for y ~ 6.0 @|||). So the 
phonon mean free path in the y — 6.0 sample unam- 
biguously reaches its maximum, boundary-limited value 
at ~ 0.15 K. It is then reasonable to expect a very sim- 
ilar phonon behavior in the y — 6.9 sample, given its 
nominally identical crystalline quality and surface qual- 
ity, and its comparable dimensions. This is nicely borne 
out by the k/T data in Fig. 1: the only difference be- 
tween the two curves (y = 6.9 and 6.0) is a rigid offset. 
In such a well-defined context, the appearence of a linear 
term upon introducing electronic carriers is conclusive ev- 
idence for the existence of zero-energy quasiparticles in 
YBa2Cu30o.9 and, as a result, it confirms a key feature 
of the basic theory of transport in unconventional super- 



conductors |)-B J14|-[rj| . In this connection, earlier claims 
of a residual electronic linear term in k of YBa2Cu307_a 
were inconclusive, being all based on the same analysis as 
used here but applied to arbitrary temperature regimes 
(for a review, see |Iql). 

TABLE II. Parameters used in fitting k/T to a + bT 2 , 
where a = kq/T is the electronic residual linear term and 
b = K p h/T [i is the asymptotic phonon T 3 term, w is the 
mean sample width and Ao is calculated from Eq. (2) using 
(3 = 0.3- 0.4 mJ/K 4 mole and (v ph ) = 4000 m/s. 
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FIG. 2. a-axis thermal conductivity of the four Zn-doped 
crystals, plotted as k/T vs T. 

Having established the existence of a residual normal 
fluid in YBa2Cu306.g, the next question is that of uni- 
versality. This is addressed by looking at concentrations 
of Zn such that T ranges from < 0.014 up to 0.54 T c . The 
thermal conductivity of YBa2(Cui_2;Zn x )306.9 is shown 
in Fig. 2, where it is apparent that k is unaffected by the 
variation in T at T ~ 0.1 K, where the heat is carried pre- 
dominently by quasiparticles (cf. Fig. 1). In other words, 
transport by the residual normal fluid is universal. 

The T ->■ limit of k/T is obtained from a fit to a+bT 2 
limited to T < 150 mK, as applied earlier, which yields 
the values for a = kq/T and b listed in Table II. Note 
that the ratio Ao/(2w/y / 7r) ~ 1 for all crystals, proving 
that the asymptotic phonon regime was reached in all 
cases. (The somewhat larger ratio for the pure sample 
is intriguing - further work is needed to elucidate this.) 
As seen from a plot of kq/T versus T, shown in Fig. 3, 
these values are consistent with a universal linear term 
of 0.19 mW K~ 2 cm -1 . Note, however, that the error 
bars on the values of a and b in Table II are fairly large, 
because they combine uncertainties on the geometric fac- 
tors (largest for the rather short 2% sample) and on the 
fit, which is limited to a small temperature range (small- 
est for the 3% sample). One way of eliminating the un- 
certainty on the geometric factor is to use the resistivity 
data obtained with the same contacts. Indeed, by fixing 
B = 1.03 fj,£l cm K^ 1 for all samples, thereby imposing 
the reasonable constraint that the inelastic scattering is 
not affected by small levels of Zn, one can correct Ko/T 
by multiplying it by B(x)/1.03. This yields the following 
corrected values: 0.17 ±0.01, 0.17 ±0.02, 0.23 ±0.02 and 
0.21 ± 0.04 mW K~ 2 cm" 1 , for x = 0, 0.6, 2 and 3%, re- 
spectively. These are plotted in the inset of Fig. 3 versus 



the similarly corrected T p . The corrected plot with its 
smaller error bars no longer allows for a constant linear 
term: there is a small but definite upward slope, with a 
minimum growth of 30% over the range of T p and a max- 
imum growth of 55%. From this we conclude that while 
the residual linear term is universal, in the sense that 
a 10- fold increase in V (from 0.014 to 0.13 T c in going 
from x = to 0.6%) leaves Kq/T unchanged, at larger T 
there is a slight increase, reaching approximately 40% at 
T/T c ~ 0.6. 
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FIG. 3. Residual linear term vs scattering rate for the four 
crystals of YBa2(Cui- a; Zn I )306.9; the dashed line indicates 
a constant at 0.19 mW K~ 2 cm -1 . Inset: same, but with 
corrected values (see text); the solid line is a least-squares fit. 

Let us now compare our results with the theory 
of heat transport in unconventional superconductors 
]i~4] [nj . The first point to emphasize is that universality 
is only expected for special gap functions with appropri- 
ate topology and symmetry. This is the case for a pairing 
state of d x 2_ y 2 symmetry, with line nodes at azimuthal 
angles 4> = tott/4 (to = 1,3,5,7). The T = limit of 
k/T along the x (or a) direction is [ ]i^ , |l5[ 



«oo 
T 



Lqo-qo 



L ne 2 2h 
m* ttS 



hkjju 2 1 
~6e^~S 



(3) 



where <7oo is the universal limit of charge conductivity 
§,!§, L a = (n 2 /3)(k B /e) 2 is the Sommerfeld value of the 
Lorenz number, and S — \dA(4>) / d(p\ noc i e is the slope of 
the gap at the node |Fq| . The topology of the excitation 
gap right at the node (e.g. A(0) ~ <jf> — 7r/4) determines 
universality and then the slope S sets the magnitude of 
kqo/T. Note that a gap function of the right topology but 
of s-wave symmetry will not in general show universal 
behavior, so that our observation of a universal kq/T is 
strong support for a gap of d-wave symmetry. 



3 



A quantitative comparison with the theory reinforces 
this conclusion. In the simplified case of the standard d- 
wave gap A o cos(20), S = 2A . Using available estimates 
of hujp and the gap maximum Ao, respectively equal to 
1.3 eV |§ and 20 meV @, one gets 



^ = 0.09 mWK^cm" 



(4) 



which is remarkably close to the measured value of 
0.19 mW K -2 cm" 1 . Given that the real gap will have 
more structure than a simple cos(2</>) dependence, the 
factor 2 discrepancy suggests that it actually rises from 
the node half as fast as in the simple model. This in no 
way detracts from the conclusion that a (generalized) d- 
wave gap is in excellent quantitative agreement with the 
universal heat conduction observed in YBa2Cu^O(,.g. 

The second point to consider in a comparison with 
the theory is the fact that universality is only achieved 
when fvy « Ao , where fvy is the bandwidth of impurity 
bound states responsible for the zero-energy excitations 
fl5|| . The bandwidth grows with T in a way which de- 
pends very strongly on the scattering phase shift 5q. It 
is largest in the limit of unitarity scattering, So = tt/2, 
where fry - y/nA hT/2 @. For the pure and 3% sam- 
ples, with Tp/T cQ = 0.014 and 0.54, this gives fry/A ~ 
0.1 and 0.6, respectively (for A ~ 20 meV = 2.5 k B T c0 ). 
Thus we expect deviations from universality for the sam- 
ples with high Zn doping. Quantitatively, the depen- 
dence of kqo/T on r was calculated by Sun and Maki 
jl4]], who find a monotonic increase, which gets to be a 
factor 1.9 at T/T cQ = 0.54 (see also Ref. 17). Such a 
large increase is incompatible with the data (see Fig. 3). 
On the other hand, a 40% growth in the residual linear 
term, consistent with the data, would agree with the cal- 
culation if r ~ r p /2, namely 0.3 T c0 for the 3% sample 
instead of 0.54 T c q. Interestingly, this is the T one de- 
duces self-consistently from the theory |^(J , based on the 
measured T c suppression. Note, however, that account- 
ing for a smaller V in terms of a smaller "effective" w 2 in 
Eq. (1) leads to an even smaller Koo/T from Eq. (3). 

These minor discrepancies notwithstanding, one of the 
main implications of the good agreement with the the- 
ory is that impurity scattering in the cuprates is well- 
described by a phase shift very close to tt/2 |p!5| , some- 
thing which has been assumed often but rarely verified. 
Nonetheless, a proper interpretation of the data should 
include the possibility of a small departure from the uni- 
tarity limit §L% . This would lower 7, making it easier to 
satisfy the condition fry « Ao, and possibly to account 
for the weak variation in kq/T. 

The present results have implications for other prop- 
erties of YBa2Cu306.9, such as specific heat JlT| , ^2] | and 
microwave conductivity pl| . It has not yet been possi- 
ble to probe the residual normal fluid reliably via these 
properties, but the upper bounds imposed by the data 
so far are consistent with the behavior predicted on the 



basis of the thermal conductivity data reported here. 

In summary, we presented the first observation of uni- 
versality in the transport properties of a superconduc- 
tor. A residual linear term kq/T = 0.19 mW K -2 cm -1 
is clearly resolved in the thermal conductivity of 
YBa 2 Cu306.g and attributed to electronic carriers. The 
observation of this residual normal fluid is a powerful 
validation of the basic theory of impurity scattering in 
unconventional superconductors. The fact that kq/T is 
universal, i.e. virtually unaffected by changes in the im- 
purity scattering rate, strongly confirms the gap has hav- 
ing d-wave symmetry. However, from the magnitude of 
Ko/T, it appears that the gap rises more slowly at the 
nodes than described by the standard function Aocos(20) 
with Ao = 2.5 ksT C Q. 
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